Anomalous strength of membranes with elastic ridges 
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We report on a simulational study of the compression and buckling of elastic ridges formed 
by joining the boundary of a flat sheet to itself. Such ridges store energy anomalously: their 
resting energy scales as the linear size of the sheet to the 1/3 power. We find that the energy 
required to buckle such a ridge is a fixed multiple of the resting energy. Thus thin sheets with 
elastic ridges such as crumpled sheets are qualitatively stronger than smoothly bent sheets. 
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INTRODUCTION AND THEORY 

When the boundaries of an elastic sheet are sufficiently 
distorted, singular deformations in the sheet often oc- 
cur H, H, ^, H, H . These singularities are immediately 
apparent when one crumples a sheet of paper between the 
hands. Analogous singularities occur whenever the sheet 
is forced to have points of high curvature j|, ^, ||] . One 
can readily form such a "puckered" sheet by fastening the 
edges of that sheet to form disclinations as shown in the 
cube of Fig. |l|. Recently it has been recognized that such 
structures typically have their energy concentrated into 
so-called stretching ridges, in which strain and curvature 
energies are m balance g |, §. The ridge increases the 
stored elastic energy in the structure dramatically rela- 
tive to that of a uniformly curved sheet occupying the 
same region. This enhancement of energy grows as a 
power of the linear size of the structures relative to the 
thickness of the sheet. 

The large energy content of these structures suggests 
that they should be strong. They should resist deforma- 
tion and buckling under external loads qualitatively bet- 
ter than nonsingular structures. In this note we demon- 
strate this strength by showing that the scaling laws gov- 
erning the ridge energy also predict its response to loads 
and its buckling threshold. 

We simulate the compression of an elastic cubic shell by 
the application of external, inward-pointing forces at its 
vertices. The cube is formed from a flat sheet as shown in 
Fig. 0. With this connectivity the sheet naturally forms 
stretching ridges along the edges of the cube. We push 
on the vertices of the cube until the ridges buckle into 
several smaller ridges. By rescaling the compressed con- 
figuration with the same thickness scaling laws as those 
appropriate for the resting configuration, we demonstrate 
that the response of the ridge to this kind of forcing is 
completely described by the 1/3 power ridge scaling so- 
lution determined in Q. We also demonstrate that the 
onset of the buckling instability follows this scaling as 
well. 

Ridge scaling comes about as a way to balance bend- 



ing and stretching energy costs. It has been shown Q 
that the important dimesionless parameter governing the 
elastic energy of a single ridge of length X and material 
thickness h is the aspect ratio 
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where Y is the 2-dimensional Young's modulus of the 
material, k is an effective bending modulus, and v is its 
Poisson ratio. If all lengths are expressed in units of 
X and all energies in units of k, the equilibrium config- 
urations of a thin sheet are completely described by a 
dimensionless form of the von Karman equations. 
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where all derivatives are taken with respect to the di- 
mensionless variables x/X and y/X. Here the bracket 
product represents 
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P is an arbitrary external pressure field acting normal to 
the surface (in units of kX^), and the potential fields are 
related to the curvature Cap and stress aap by 



Cap — X ^dadpf 

= KX^'^Cap.ep^d^d^x- 
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This form of aap automatically satisfies the equilibrium 
condition daC^ap — 0. 

Since A comes into the von Karman equations multiply- 
ing the stess source term, the possible configurations of a 
thin elastic sheet are well described by a stress free, A = 
folding solution plus boundary layers at the fold lines. 
Lobkovsky's insight in |Q was to try a scaling solution 
for the boundary layer of a single ridge which matched 
the / scaling of the outer, sharp fold solution. For a fold 
of dihedral angle a across the line y = 0, f = a \y/X\. 
Accordingly, on the boundary layer / should scale with 
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the same power of A as the dimensionless transverse coor- 
dinate y/ X. Substituting a scahng form into the dimen- 
sionless von Karman equations with P = and equating 
the highest order terms yields scaling of the form 



/ = A^/^; X = A-'/^x, V = X'^'Xy, x = X 



X, 



(6) 



where the tildes denote dimensionless, scale free coor- 
dinates and functions. This translates to A^/^ scaling of 
the boundary layer width, A~^/^ scaling of the transverse 
ridge curvature, and A^/^ scaling of the strain along the 
ridge length. 

External forcing applied to the sheet enters the von 
Karman equations via the term P and via boundary 
conditions at the sheet's edges. In this research we con- 
sider an external potential which essentially applies point 
forces to either end of a ridge. Since the spatial extent of 
the applied force is a delta function we do not expect it to 
destroy the spatial scaling of the ridge solution. There- 
fore we may reasonably expect to find that the equilib- 
rium configuration of a ridge under a given compressive 
force is identical to rescaled configurations of ridges with 
different material thickness and properly rescaled exter- 
nal force magnitudes. 

To calculate the proper rescaling of the forces on the 
vertices for a similarity solution, we consider our forcing 
as a boundary condition consisting mainly of an in-plane 
point force. This force amounts to a point stress at the 
edge of the sheet with the form 
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So, to find similar scaled configurations of the sheet, we 
must scale the same way axx scales on the ridge. 
Since j^x scales as A^/^, and 

then axx ~ f^^^^^^X^'^^xx, where "/xx is a dimensionless, 
scale free function. To express ai°J in similar fashion, we 
substitute the scale free j/- variable y = X~-^/^y/X , so that 
S{y) = \~^^^X~^S{y). Thus, the proper scale free force 
can be written in terms of Fo as 
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For reasons which will become clear in the next sec- 
tion, we cannot measure the force we apply to our ridge 
with very good accuracy. However, we can measure the 
inward displacement A of the ridge ends caused by this 
forcing. These two quantities may be related by assuming 
that the work done by equivalent rescaled forces, given 
approximately by FqA, scales the same as the total en- 
ergy of the resting ridge configuration. The total energy 
of a resting ridge scales as kX~^/^, so equivalent values 
of FqA/k will scale as A^^/"^. Given the scaling of Fo 
from Eq.||, the scale free A must be related to the actual 
displacement by 
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FIG. 1: Typical elastic sheet used in this study. The top 
image illustrates how a flat sheet is joined at its edges to form 
a cube. The bottom image shows the resting configuration of 
the cube with no external forces acting when the thickness 
of the sheet is .0004 of the the edge length and the Poisson 
ratio is 1/3. Darker shading represents higher strain energy 
density. The energy was minimized in the outlined diamond- 
shaped region; the rest of the cube shape was inferred by 
symmetry. Slight numerical symmetry breaking between the 
left and right sides of the diamond created slight mismatches 
of the inferred surfaces on other faces, such as the right-hand 
face. The numerical grid is visible as a quiltlike texture. It 
has a finer scale at the edges and corners where curvature is 
larger. 



NUMERICS 

To test these scaling predictions we simulate a fiat 
elastic sheet joined to form a cube as shown in Fig. [|. 
The sheet is a triangular grid with variable grid spacing. 
Strains and curvature are taken to be constant across the 
face of each triangle, with curvature calculated relative to 
the local normal at each triangle. Bending and stretch- 
ing energies were assigned to the curvature and strains on 
each triangle using the forms for elastic energy presented 
in j^] for a sheet of elastic thickness h and Poisson ratio 



A^/^xA. 
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FIG. 2: Energy of ridges at rest and at the buckling threshold. 
Straight lines are least squares fits to a scaling form y — ax''. 
In this plot A ranges from 1.25 x 10"^ to 1.77 x 10"^ The plot 
shows the total elastic energy {Eb + Es) in the sheet after 
minimization. The scaling exponent fit for the resting ridge 
values (lower line) was —0.32, the fit at the buckling threshold 
was —0.31. The inset shows the difference between threshold 
energy and resting energy in units of the resting energy. This 
energy ratio is best fit by a scaling exponent of 0.05 ± 0.02 
and is consistent with a constant ratio. 

1/3. The gridding was chosen to have smooth gradients 
in triangle density over most of the surface while concen- 
trating the lattice spacing at the vertices by a factor of 
10'^ and across the ridge-line by a factor of 10^ compared 
to the flat regions far from the ridge. The concentrations 
factors were chosen arbitrarily, within the limits of the 
mapping, to make the gridding near the vertices as fine 
as possible. 

Pushing on the tips of the cube is accomplished by 
introducing repulsive potentials of the form V{r) — 
Cp/ Ir — Xp\^ centered around points Xp which lie just 
outside the vertices of the cube. To prevent rotation, the 
vertices are explicitly constrained to lie on radial lines 
that pass through the center of the cube and the points 
(±1, ±1, ±1) which define the vertices of a perfect cube. 
The center points of the pushing potentials are located 
on these radial lines, at the point where the vertex would 
lie if the sheet were sharply folded at the ridges - re- 
laxation of the ridge curvature draws the vertices inward 
from these points for an unforced resting ridge. This 
potential concentrates the forces at the corners without 
creating lattice-scale numerical instabilities. Cp was var- 
ied to apply different loading. 

An inverse gradient routine |^ was used to minimize 
the total elastic and potential energy of the sheet as a 
function of the coordinates of all the lattice points for 
given parameters k,Y and Cp. 

To save computational time, energy was computed and 
minimized on only one ridge- line (the diamond shaped re- 
gion highlighted in Fig. |l|) . The positions of points on the 
rest of the cube were calculated by reflection across sym- 
metry planes. This constraint explicitly required that 
the edges of the simulated region be confined to the sym- 



metry reflection planes. Thus many forms of deformation 
that break the symmetry of the resting state are not pos- 
sible in the simulation. However, the minimization may 
break the left-right symmetry of the diamond-shaped re- 
gion, and we found that it does so to a slight degree. The 
assymmetry shows up in our reconstruction as a slight 
mismatch in other faces of the cube, as noted in Fig. 
We do not believe these mismatches are important for 
the scaling phenomena we report. 

We found minimum energy configurations for ridges of 
aspect ratio A ranging from 1.25 x 10"'^ to 1.77 x 10~^. 
The upper bound on A was determined by the range of 
validity of the ridge scaling solution - above this value 
the width of the ridge becomes comparable to that the 
sheet. At the other extreme, for A < 10^^ the radius of 
curvature at the ridge line becomes comparable to the 
spacing of our lattice and the simulation ceases to be 
accurate. For each value of A we first found the mini- 
mum energy configuration with no applied forces. This 
configuration was then used as the initial condition to 
find minimum total energy configurations for the ridge 
in the presence of the tip-pushing potential described in 
the numerics section. We made the tip-forcing progres- 
sively stronger by increasing the potential coefficient Cp 
in constant steps, each time evolving the grid to an en- 
ergy minimizing configuration. The forcing was increased 
until the ridge buckled. Buckling of the ridge was marked 
by the appearance of new points of sharp curvature along 
the ridge as well as a sudden decrease in the total elastic 
energy Es + Eb of the sheet. Thus the unbuckled sheet is 
metastable. We made no attempt to scale the size of the 
potential step with the thickness of the sheet, so while it 
took nearly 30 steps to buckle the ridge in thicker sheets, 
it only took 6 or 8 steps for the thinnest sheets. 

FINDINGS 

The plot in Fig. ^ shows scaling of the total elastic 
energy in the cube versus A for ridges at rest {Cp — 0) 
and at the buckling threshold. Scaling of the total elastic 
energy for the resting configuration is consistent with a 
kA~^/^ dependence, in agreement with prior theory and 
simulation ^ . Fig. |^ also shows that the elastic energy 
measured at the buckling threshhold exhibits exactly the 
same scaling as on resting ridges. This suggests that 
the response of the cube to our tip-forcing is completely 
determined by the ridge, and that the particular form of 
our forcing potential docs not destroy the length scaling 
of the ridge. The inset in Fig. || shows that that the 
energy correction at the buckling threshold is nearly a 
constant fraction of the total ridge elastic energy. 

Scaling of the force response is verified by the demon- 
stration presented in Fig. ^ of a similarity solution for 
ridge shape. Plot (a) in this figure shows values of Cyy, 
along a line in the material coordinates which bisects the 
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FIG. 3: Demonstration of a similarity solution for the ridge re- 
sponse to forcing. Both plots show Cyy, the curvature across 
the ridge-line, versus the y material coordinate on the line 
which bisects the ridge-line. The data is for sheets with seven 
different values of A, ranging from 1.25 x 10"^ to 1.25 x 10"". 
Plot (a) shows Cyy X (A/Ai)^''^ vs. y x (A/Ai)"^/^ for the 
ridges at rest and for ridges with inward vertex displacement 
A(A) = A^"' x (A/Ai)"'^'', where A is measured from the 
resting vertex positions, Ai is the aspect ratio for the thick- 
est sheet, and A^^'^' was vertex displacement at the buckling 
threshold for the thickest sheet. The profiles with the large 
central peak are the buckling threshold values. (The small 
dimple in the data at y = is a numerical artifact due to 
a discontinuity in the gridding density across the ridge-line. 
For finer gridding this dimple goes away, while all other local 
values of curvature remain constant.) Plot (b) shows unsealed 
Cyy versus y for the buckling threshold profiles plotted in (a) . 



simulated ridge line, for several different sheet thicknesses 
and two different values of rescaled ridge tip displacement 
A. We tested values of A which spanned the ridge re- 



sponse from zero forcing to the point of buckling. The 

1/3 

magnitude of Cyy was rescaled by (A/Ai) ' and the y 
coordinate was rescaled by (A/Ai) ' . We found that 
the rescaled Cyy profiles mapped onto each other best if 
the values of A used for each sheet thickness were de- 
termine by rescaling Ai, the value used for the thickest 
sheet, by (A/Ai)"'^^, where Ai is the aspect ratio of the 
thickest sheet. For comparison, the unsealed Cyy versus 



y for a particular rescaled A is shown in Fig. ^(b). The 
rescaling exponents used for Cyy and y were taken from 
the ridge scaling solution. The A rescaling exponent is 
very close to the theoretical value of 2/3 derived above. 

We have demonstrated scaling of the ridge response 
to a typical form of external forcing. Since the buckling 
threshold of the ridge under external forcing is deter- 
mined by the energy scaling of the resting ridge solution, 
the relatively high elastic energy content of very thin 
ridges translates directly into increased ridge strength. 
We are actively pursuing a better understanding of the 
buckling properties of this distinctive elastic structure. 
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